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Abstract 

A Fortuin-Kasteleyn cluster on a torus is said to be of type (a, b), a, b e Z, if it possible to draw a curve 
belonging to the cluster that winds a times around the first cycle of the torus as it winds — b times aroimd 
the second. Even though the Q-Potts models make sense only for Q integers, they can be included into a 
family of models parametrized by |3 = for which the Fortuin-Kasteleyn clusters can be defined for any 
real |3 € (0,2]. For this family, we study the probability 7r({a, b}) of a given type of clusters as a function 
of the torus modular parameter t = Tt -I- ixi. We compute the asymptotic behavior of some of these 
probabilities as the torus becomes infinitely thin. For example, the behavior of 7t({l , 0}) is studied along 
the line Tr = and Ti — > oo. Exponents describing these behaviors are defined and related to weights hr,s 
of the extended Kac table for r, s integers, but also half-integers. Numerical simulations are also presented. 
Possible relationship with recent works and conformal loop ensembles is discussed. 
Keywords: Fortuin-Kasteleyn clusters, torus, homology probabilities, homotopy probabilities, percolation, 
Ising model, logarithmic minimal models, SLE. 
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1 Introduction 



One of the main observables of two-dimensional percolation is the crossing probability between two dis- 
joint subsets of the boundary of a domain. This domain is usually taken homeomorphic to a disk. As Lang- 
lands and his colleagues [11 J were finishing their numerical study of universality and conformal rnvariance 
of crossing probabilities, 1. Gelfand suggested to explore percolation on compact Riemann surfaces. The 
simplest surface to study is the torus and the most natural observable is then the homologic properties of 
the percolating cluster, or more precisely, the probability that a configuration contains a homologically non- 
trivial cluster. (Since these clusters are geometric objects, it might be easier to think about their homotopic 
properties instead of their homological ones.) Let cui and 0)2 be the two-dimensional linearly independent 
vectors along the two sides of the parallelogram defining the torus. In the following these will be identified 
to points in the complex plane. If a non-trivial cluster exists and if it winds a times along cui of the torus 
while it wraps b times along — cui, the cluster is said to be of type {a, b}. All other non-trivial clusters of that 
configuration, if any, will be of the same type. (The integers a and b are coprimes. Types {a, b} and {—a, — b} 
are considered identical.) For that reason, the homology property of a configuration may be defined as the 
type of its non-trivial clusters. If the configuration contains no non-trivial cluster, it is said to be of type {0}. 
Finally, if the configuration contains a cluster that has both a path around the first cycle, that is along cu 1 , 
and a path along CU2, this configuration is of type Z x Z. With that notation, each configuration is associated 
with one of the subgroups H of the homology group Z x Z of the torus: {0}, Z x Z and {a, b} with a, b co- 
primes. The same notation {a, b} is used for the type of a configuration and the subgroup generated by an 
element of that type. Langlands et al measured the probability of a few of these subgroups for percolation 
and gave some numerical evidence for their conformal invariance. 

Pinson [17] obtained analytic expressions for the probability of these various subgroups as functions 
of the quotient t of the fundamental periods cui , W2 G C of the torus. His computation relies on a clever 
argument giving an orientation to the curves bounding clusters. (See IIT4l l9l.) This is done in a way that 
does not change the partition function, but does allow for the identification of the homology properties of 
intervening clusters. His computation is mathematically rigorous, except for the step taking the limit as the 
mesh goes to zero; for this, he used Nienhuis' renormalization group argument [14J that ties the quantities 
under study to known results for the Coulomb gas. A more rigorous treatment of this step remains open. 

Arguin [IJ extended Pinson's argument to Q-Potts models, Q — 1,2,3,4. To do so, he considered the 
Fortuin-Kastele3m graphs or clusters of configurations. These are the natural extensions of the clusters of 
percolation, the Potts model with Q = 1 . Arguin showed that Pinson's formulae need only a small change 
for the Q-Potts model with Q > 2. He also supported his new expression with numerical data for the four 
integer values of Q. 

Works on or using probabilities of homology subgroups of FK clusters has not been limited to the the- 
oretical predictions. Ziff, Lorenz, Kleban [23| were the first to provide numerical support for their univer- 
sality. Later Newman and Ziff [13 1 used them to give a precise estimate of the critical probability for site 
percolation on a square lattice. It was then the most precise available estimate. And recently they were 
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again used to obtain precise estimates for critical probability for percolation on several lattices |I7|. (These 
probabilities are called wrapping probabilities in these works.) 

In the definition of Potts models, Q gives the number of states accessible to the basic random variables, 
often called spins. As such, Q must be an integer When the partition function is rewritten in terms of 
Forturn-Kasteleyn graphs (hereafter FK graphs), the parameter Q appears in the Boltzmann weight as Q'^'^ 
where N c is the number of FK connected components in the configuration. In this formulation, the condi- 
tion that Q be an integer may be relaxed. One then gets a one-parameter family of models, usually studied 
for the values of Q in the interval (0,4]. It is between this family of models and the family of stochastic 
Loewner processes that a close tie seems to exist, and has been established for some particular cases. The 
stochastic Loewner equation with parameter k (SLEk) is believed to describe the growth of the boundary 
of a FK graph. The exact relationship between the two parameters Q and k is 



with KG [4, 8) and, again, Q e (0,4]. Percolation corresponds to k = 6 (and Q = 1) and the Ising model to 
^ — ^ (Q — Thfi mathematical tools to describe not only the boundary of a single FK cluster, but the 
set of loops described by the boundary of all clusters in a configuration are now emerging. Conformal loop 
ensembles, defined by Camia and Newman for percolation [4J and more generally by Werner [22J (see also 
Il3l). might allow for the rigorous study of homological properties of configurations, as defined and studied 
by Langlands et al, Prnson and Argurn. 

The goal of the present paper is to extract from the known expressions of the probabilities for the vari- 
ous homology subgroups their asymptotic behavior for two limiting cases. The first is when the quotient t 
of the periods goes to infinity or to a real rational number The second is when Q goes to zero. The reason 
to study the latter is mostly curiosity. For the former, the reason is twofold. Many results proved using SLE 
techniques describe asymptotic behavior The first reason is therefore to seek exponents to describe limiting 
behavior that might be easier to obtain with SLE (or CLE). The second reason is to probe deeper the rela- 
tionship between SLE and conformal field theory (CFT). Several critical exponents appearing (rigorously) 
in the context of SLE had been predicted within CFT, and a large subset of these appeared in the Kac table 
of the associated minimal conformal model. It is agreed, but not proved, that SLEk describes properties of 
the conformal theory with central charge 



Minimal models appear when c and k are rational. Let k be rational and of the form 4p '/p with p ' > p > 1 , 
coprune integers. The conformal spectrum of the minimal model with central charge c = c( k) is constructed 
from the Virasoro highest weights 



It has been recognized however that the minimal models, constructed out of finite sets of primary fields and 
therefore of highest weights Hr , s , are probably too restrictive and might not capture all physical observables . 
Half-integers r and s have been considered [20J and several works about logarithmic minimal models have 
shown that the upper bounds on r and s need to be relaxed. (See, for example, lil2i,15| for recent arguments.) 
Maybe one of the most striking examples of this fact is Cardy's formula that describes the probability of 
crossing within a rectangle for percolation. For limiting geometries, that is for rectangles very wide or 
narrow, the probabilities approach or 1 with the power of h] 3 = 5/ ari exponent that does not belong to 
the minimal set. Another example is related to the problem studied in the present note. In |]2j, Arguin and 
Saint- Aubin showed that, when the quotient t of the fundamental periods of the torus tends to zero along 
the imaginary axis, the probability 7t({1 , 0}) for the Ising model goes to 1 as intuitively it should, but more 
precisely it goes as 7t({l ,0}) — > 1 — (q^) sfi (q'^) — (q^)Tf2(q^) — . . . where q — e^"'^ and f 1 and fi are analytic 
in a neighborhood of q = 0. The exponents are twice the highest weights H] 2 = yg and hs 3 — ^; the first 
belongs to the spectrum of the minimal model, the second does not. It is this observation that led us to ask 




-r,s — 



(Kr-4s)2 - (k-4)2 
16k 



l<r<p-1, 1<s<p'-1. 
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Figure 1: The torus in the complex plane, with t = 



-2/3 + 1 



whether exponents obtained by taking limits of the geometry would always be in the extended Kac table of 
the corresponding models when k is rational. (Every conformal weight hr s is repeated an infinite number 
of times in the extended Kac table. Arguin and Saint-Aubin chose (r, s) — (1 , 2) and (3, 3) for the leading 
exponents of the Ising model. We shall come back to this choice after determining the exponents for the 
general case.) 

Our notations are the following. The torus is identified with the quotient C/jcui , cua) where {tui , wz} is 
the integral lattice generated by cui , wz G C such that 0, cui , cu2 are not colinear. We choose cui = 1 and 
Im CU2 > 0. Their quotient t — cui/cui is the modulus of the torus with Tr and Tt > its real and imaginary 
parts. Figure [l] shows these basic elements for a torus with t = — | + i. We follow the convention set in 
IITTI ni for the winding numbers: they are positive in the direction of cui and — CU2- Figure |2] shows FK 
configurations of three different types drawn on the torus T = i. Configuration (c), for example, is of type 
{2, — 1 } according to the above convention. 

It is natural to break the partition function into sums over configurations of a given tj'pe or generating 
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Figure 2: Examples of FK configurations of type (a) {0}, (b) Z x Z and (c) {2, — 1 } groups, drawn on the torus 
with T = i 



a given subgroup H. If a A b denotes the greatest common divisor of a and b (with a A = a for all a), the 
partition function is 

Z = Z({0})+Z(ZxZ)+ Y_ Z({a,b}). (2) 

aAb = l 

The observables under study are the probability of a given subgroup H, namely 7t(H) — ^ . All these 
quantities depend on the size of the lattice covering the torus and the model labelled by Q. (For clarity we 
sometimes add an index, Q or t, to quantities under study, e.g. Z — Zq .) Their thermodynamic limit, when 
the mesh size goes to zero, are known at the critical temperature. The expressions obtained by Pinson IflTl 
for Q — 1 and generalized by Arguin 1 1 1 for Q e {1,2,3,4} are 



where 



Z({Q,b}) = ^ Zbk,ak{g/4)(cos[7teolc] - cos[7tk]) (3) 

keZ 

Z({0}) = 1 Y. Z„,,^,{g/4)cos[7t(mAm')] (4) 

m,m''eZ 

Z(Z X Z) = Q X Z({0}) (5) 

(6) 



and 

Q =4cos2[7Teo/2], g=4-2eo, eo=2-8/K, Q G (0,4], eoe[0,1), kg [4,8). (7) 

The parameters Q , g , eo and k are all in one-to-one correspondence to one another in their respective range. 
(We use them in the way historical developments have introduced them.) Dedekind function is r|(q) = 
qi /24 — q'^). Pinson's and Arguin's arguments extend trivially to the models of Fortuin-Kasteleyn 

cluster with a real Q in the interval (0,4]. We use these expressions as our starting point. 

The paper is organized as follows. In the next three sections, we study the following three limits: of 
7t({1 , 0}) when t = ixt and Xi — > oo, of 7T({a, b}) when i — -g + iTi, with Xi —> and finally of 7t(H) for any 
H c Z X Z when Q — > 0. The last section is devoted to Monte Carlo verifications of some of the results. 



2 The probability 7r({l , 0}) in the limit = 0, Ti ^ oo 

The first limit to be studied is when t = ixt with Ti — > oo, i.e. the limit when the torus becomes a very 
thin ring. The corresponding parallelogram in the complex plane becomes an infinitely tall rectangle of 
constant width equal to 1 . Curves winding once along w i become very likely. In fact their relative length 
with respect to those winding once in the direction uj2 becomes negligible and it is therefore expected 
that, in this limit, all configurations will have curves of t5rpe {1 , 0} and none of type {0, 1}. In other words, 
7t({1 ,0}) — > 1 and the probability of all other groups goes to 0. What should be the expected behavior of 
7t({l,0}) for finite but very large ii? Cardy's formula [51 provides a fair guess. This formula gives, for 
percolation, the probability 7th of horizontal crossing in a rectangle of width H and height V as a function 
of the aspect ratio r = V/H. For luniting geometries the probability behaves as 

7th{r) — > cie-"/^'" and 1 -7TK(r) — ^ cie""''/^ 

T— iO r— >oo 

for known constants Ci and cz- Even though the intersections of a percolating cluster with the left and 
right edges of the rectangle might be in general at different height, these two intersections are likely to 
have points with the same vertical coordinates if the rectangle is very narrow, that is when r ^ 0. Such a 
percolating cluster would be a FK cluster of type {1,0}, if opposite edges of the rectangle would be glued 
together. Therefore one may expect the following behavior 

7t({1,0}) = 1-^c^qT" (8) 

n 
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with positive exponents yn and the natural parameter q — e^^""^' if the real part of t vanishes. Note that q 
goes to when ti — > oo. The goal of this section is to determine the leading exponents yn as a function of Q 
or, equivalently, eo- (Some care should be exercised as the immediate extension of q to a t in the upper-half 
plane by q = e^"^'^ does not coincide with the usual definition of the nome of elliptic fimctions which is 

The probability 7t({l , 0}) is given in the form Zq ({1 , 0))/Zq . The first step is to express the numerator 
and denominator in a form suitable to extract these exponents. From ||3|: 

Zq({1,0}) = y Zok,ik{g/4)(cos[7teolc] -costTtk]) = , \,y \l^y e""^ (cos[7Teok] - cos[7Tk]). 
— m q V 4Ti ^ — 

1 

To rewrite the e in terms of q, Poisson summation formula will be necessary: 

neZ ^ keZ 



After expanding the cosines in terms of exponentials, Poisson formula gives 

1 



^ " keZ 



Since the function q^^/^'*ri(q) has a Taylor expansion, the above form allows for the identification of the 
leading terms in the numerator. Note however that the expansion of |r|(q)P will not be used, since this 
same factor appears in the denominator. 
The denominator 

Zq = (Q + 1)Z({0})+ Y_ Z({a,b}) 

aAb = l 

has two parts, which will be tackled separately. The partition function restricted to configurations with 
only trivial clusters is 

^'3«0}) = 2hiW\/4 ^ e-'^'^^i^cosWmAm')]. 
To get rid of the cos[7T{Ta A m')], we notice that 

Z = Z +( Z - Z ) (11) 

m,rrL'GZ m,m'e2Z m,m'£Z m,m'G2Z 

In the first sum, both ra and m' are even which makes ra A ra' even and cos[7T(m. A ra')] — 1 . The other 
terms, in the parenthesis, are terms for which either m or ra' is odd, and cos[7T(ra A ra')] = — 1 . Therefore: 



zq({0}) = — ^./^(2 y y_ 

4riqKVTi\ ^ — ^ — 

V L ^ Tn'c7Z. Tn m'f: 



rri,m'G2Z m,m'e2 

Sums over multiples of an integer f e N will appear often and it is useful to define 

c^(f,g)-\/^ Z 

V ^ — 

* m.m'efZ 



Q r _ 7T9_f£_m/^ Tigf^m^T^ 

y I \ „ a.. \ g , 



- Y_ q^ + ^ (12) 



where Poisson formula ||9| was used again in the last line. The partition function Zq ({0}) is then 



The remaining term of Zq, that includes configurations with non-trivial FK clusters of type {a, b} for all 
a and b coprimes, is more complicated. The sum 

Y_ ZQ({Q,b})= Y_ Zm,m'(g/4)(cos[7Teo(mAm')] -cos[7t(mAm')]) (14) 

aAb = l m.-tn'eZ 

contains two terms. The second with cos[7T(ra A m')] is exactly twice the partition function Zq({0}) just 
calculated. The first with cos[7Teo(m A ra')] does not simplify as easily; the sums must be reorganized 
before l|9) is used. To do so, consider, for m fixed, the function cos[7Teo(m Am')]. When m is non-zero, it is 
periodic in m' with period ra. Therefore 

^ Z^,^.(g/4)cos[7teo(mAm')] =^ Y_ C(d (g/4), m^O (15) 

m'eZ dlmm'edZ 

with 

C(d,eo) = Y_ cos(d27Teo)H( — ) (16) 

where ii.(x) is the Mobius function of x. (Recall that — 1 , |j.(n) = if n has repeated prime factors and 
|j.(n) — (—1 Y if TL is the product of I distinct primes.) To get | 15p6 1, the sum over m' was divided into sums 
over subsets which have the same value of cos[7Teo (ra A ra')]7in a fashion similar to the splitting proposed 
in equation (TT) . These subsets are closely related to the divisors of m, therefore leading to the splitting into 
sums over the multiples of these divisors. We must stress, however, that the only divisors to be considered 
in d|ra are the positive ones. The remaining sum can be written with the help of (TSll as 



^ Zm,m'(g/4)cos[7teo{mAra')] = ^ Zo,m'(g/4) cos(7Teom') ^ ^ ^ C{d, eo)Zn^,m'(g/4) 

m.m'GZ m'eZ mGZ* d|m m'GdZ 

= Y_ Zo,m'(g/4)cos(7Teom') ^ ^ ^ C(d, eo)Zm,dm'(g/4) 

m'GZ mGZ* d|mm'eZ 

where Z* = Z \ {0}. In the above expression, the terms with m = get a special treatment because of 
the particular definition of m A ra' when ra is 0. These were already encountered in the computation of 
Zq {{1 , 0}) and are equal to 

Y_ Zo.nx'(g/4)cos(7reom') = — i-^^q^C^+^o/iiVg. 

For the triple sum, the sum over divisors can be rearranged using 

Y_ ^h(ra,d)= Y_ Y. Mm, d) = Y_ Y. Mmd, d) 

meZ* dim deN' medZ* meZ* deN* 

and similarly for the sum of di | d in C ( d, eo ) . These manipulations have doubled the number of sums in (Ti) 
from two, on ra and ra', to four, on ra, ra', d, 6.2- This is the price to pay to use Poisson formula on the sum 
over ra' and cast everything into powers of q. The result is 

Y_ Zm,„^.(g/4)cos[7Teo(raATa')] 

= X Z Z C{d,eo)Zdm,dm'(g/4) (17) 

meZ* deN' m'ez 

1 ^ ^ ^ cost7reod2)H(d) ^^l^ + ^ ^^^^ 



and the complete partition function Zq is 



,'2 . „m2 



keZ m,m'eZ 



+ y cost7reod2)^(d) ^^— ^ + 



mez* 

d,d26N* 
m'eZ 



The probability 7t({l , 0}] is the quotient of Zq ({1 , 0}) given in |[T0) and of Zq . 

It is now straightforward to see that the lowest-order term in q is |^ for both the denominator Zq and 

the numerator Zq ({1 , 0}). After simplification of the common factor q'^f^^^, an expansion can be done to 
obtain the whole sets of exponents. An exhaustive list of possible exponents is given by taking exponents 
in the numerator and in the denominator, plus any integral linear combinations of them which arise from 
higher order terms in the expansion. The possibility that some of them could have vanishing coefficients is 
not excluded. 

It is interesting to compare the leading exponents with values Q given by CFT in the Kac table I|8|- In 
terms of g and Cq they are 

K,. . (20) 

for r, s positive integers. Note that |^ is half the power of q that was substracted to simplify the numerator 
and denominator. The first exponents for 7t({1 , 0}) are given by 

Yi = — f-r> Y2 = ^ — — (21 
4(2- eo) 2-eo 

and their integer multiples. On the range of eo, yi > Ti ■ The two exponents become equal in the limit 
eo = 1 (Q = 0); this particular case will be studied in section|4] 

Coincidences of these leading exponents or higher ones with elements from the Kac table, if any, will 
occur in the form y — 2hr,s for some r, s because of the contribution of holomorphic and anti-holomorphic 
sectors. Such coincidences do occur. The simplest r and s giving yi are r — 1, s =0 and, those giving 
y2, r = 0, s = 1. It is somewhat unusual to choose vanishing s or r. Recall however that, for logarithmic 
minimal models, the Kac table is extended and the periodicity of elements hr,s — ^r+p,s+v' model 
with K = 4p'/p allows to choose r and s positive. For some minimal models, it is however impossible to 
account for yi with integers r and s. Half-integers must be used. Arguin and Saint-Aubin [2J identified the 
two leading exponents for the Ising model to 2h] 2 and 2h3 3. Note that, when either r or s is zero, then 

s = _s. Moreover, if half -integer indices are included, the periodicity property can be refined to 
h-T.s = h^+p/2,s+p' /2- The Ising model corresponds to p = 3,p' —4 and their exponents are related to ours 
by Hi ,2 =h_i/2,o = hi/2,0 and h3,3 = ho -i =Ho,i. 

These two exponents yi and y2 are related to the fractal dimensions of geometric objects, namely the 
mass and the hull of a cluster respectively. (See [21 , 10 20J. For an extension of these geometric objects to 
loop gas models, see |19|.) In the FK formulation of the Q-Potts models, the FK cluster mass attached to a 
site is the number of bonds in the component of the FK graph containing this site. In the plane, the hull of a 
FK cluster is the set of bonds that can be reached from infinity without crossing any bond from the cluster. 
(On a torus, each cluster has an inner and an outer huU.) Their fractal dimension is 2 — 2A where A is hi /2,o 
for the cluster mass and Ko j for the hull. 

A natural explanation for hi /2,o in the present context is provided by Cardy IH (see also II20I ). Note 
first that the only way to keep a configuration from having a cluster of type { 1 , 0} is to have a cluster in the 
vertical direction. It is likely that its type will be {ra, 1 } for some ra G Z or Z x Z. Cardy gives an expression 
for the probability P{n, k) of having n clusters connecting the two extremities of a cylinder whose length is 

k times the perimeter of its section. He finds log P(n, k) 3^ ^ ^)k if n > 2. He points out that this 

expression evaluated at n = 1 is not the probability of having a single cluster between the two extremities. 



but it is the probability of having a single cluster between the extremities that does not wind in the other 
direction. When all configurations with a single cluster are considered, disregarding their behavior in the 

other direction, the probability is larger and given by log P(l , Ic) |jk. Because Zo — j for percolation, 

our first correction term is q^' = g-57rri/24^ agreement with his result. 

In a recent study of percolation, Ridout [18 1 has argued that the primary field responsible for changing 
boundary conditions in the computation of Watts' formula should be 4^2,5/2- This identification forces him 
to shift, in the extended Kac table, the admissible values of s by j when r is even. One would like to see a 
relationship with our identification of yi as 2hi /2,o- However it is r that takes an half-integer value in our 
case, and s in his case. Moreover the value hj /2,o — ^ does not appear in his shifted extended Kac table. 

The other exponents in the numerator of 7t({1 , 0}) are also part of the extended Kac table. They appear 

with r = k + 1 , s = 1 for and r = k + 1/2, s = for Not all the exponents 

of the denominator however appear in the extended Kac table, even if one allows half -integers r or s. For 
example the denominator of the exponents ^ ^2 appearing in the last sum of Zq is not bounded. There 
is no hope to find them all in the extended Kac table. Could these terms drop out of the sum because of 
cancellations? The general case is difficult to assess, but this happens in simple cases. 

It is indeed possible to find simpler form for the denominator for the four integral values Q — 1 , 2, 3, 4 
that correspond to eo — f , j , 5 , 0- For d > and these particular values of eo and for eo = 1 , the function 
C(d, eo) is particularly simple: 

C(d,0)=6d,i (22) 
C(d, ^) = - 6d.2 - ^ + 36d,6 (23) 
C(d,^) =25d,4-6d,2 (24) 
C(d, -) = ^^ 2^ 

C(d,l)=25d,2-6d,,. (26) 

The proof of these formulae is given in the Appendix. The contribution of Y.mm'ei, ^Tn,m'{g/4) cos[7teo{TaA 
m')] to Zq is then much simpler. It is 



-^(q"#+3-'-q^ + ^) forQ = l, 1 ^ (q^+e-^ - q^ + ^) for Q 
^^(q^ + i5„^^_q^ + i^_q^ + i^^qi^ + ^j for Q = 3 

and ^ q 2^™ for Q = 4. 

All the sums above are on ra, m' € Z. It is then straigthforward to show that, upon simplification of the 
factor q'^o/^9, these forms (and therefore Zq) can be written as a product j (q^^'fi(q^))(q^^' fj (q^)) 
of two finite sums where ft, fj are analytic in a neighborhood of and where all the hi, hj belong to the 
corresponding extended Kac table for some r and s integers in the range [0,p]. 



2 
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3 Probabilities 7r({a, b}) in the limit Ty = c/d, Ti ^ 

The probabilities 7t({a, b}) — ({a, b}) also have a limit, either or 1 , when t approaches a rational number 
on the real line. To see this first intuitively, consider the case {a.b} — {—2, 1}. Figure [3] presents two tori 
whose modulus parameter is on the line —2 + iTi. For each, two neighboring fundamental parallelograms 
have been drawn. A curve linking the origin to the vertex atz — Ivxi, like those shown, is of type {—2, 1 }. A 
curve of type {—2, 1 } does not need to start at a vertex, of course, but those drawn show how the curves of 
this type will be prevailing. Indeed, as the modulus parameter t slides down the vertical line = —2, these 
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Figure 3: Curves of type {—2, 1 } become more likely as t approaches t 2 + iO 



curves become very short and likely. Therefore the probability 7Tt:({— 2, 1}) should converge to a number 
larger than when t — > —2. This section shows that it actually goes to 1 . 

We have identified a torus with its modulus t, a complex number in the upper-half plane H. As it is 
well-known, this correspondence is not unique, since any pair cuj — racu] + ncua and iu^ — pcui -|- qcui 
with m, n, p, q e Z and raq — np — 1 describes the same torus, but with a new modulus t' = 0)2/0;^' . The 
special linear transformations { ^ ^ ) with integer coefficients and determinant 1 form the modular group 
SL(2, Z). It is generated by two matrices 

j ^t) ^'^ *=(o 1 

whose action on t is 

tA-1/t and tiAt-M. (27) 

The probabilities 7t({0}) = 7Tt:({0}) and 7t(Z x Z) = Tt-rfZ x Z) are invariant under the change of t by an 
element of SL(2, Z), but the probabilities 7TT:({a, b}) are not. Arguin fl] gave their transformation laws 

7T^({a,b}) =7t^+,({a + b,b}) =7t_i/T{{-b,a}) (28) 

or, equivalently 

7r,({Q,b}) =7Tg,(g.{Q,b}), geSL{2,Z) (29) 

where t h- > gx denotes the action defined by | [27| and g • {a, b} stands for the matrix multiplication g ( b )■ 
These transformations follow immediately from the form ||3} of the partition function Zq ({a, b}). 

A simple application of the modular transformation gives 7t({0, 1}) in terms of 7t({1,0}), namely 7Tt({0, 1}) = 
7t-i/T({1 , 0)). The result of the previous section implies easily that 

7r,=o+i., ({0, 1 }) = 7T,=o+i/., {{1 , 0}) = ({1 . 0}) 



where the partition functions are evaluated at q' = e^^"^^'. The limiting behavior X\ — > 0+ will therefore 
be characterized by the same exponents obtained for 7t({1 , 0}) when — > 00. 

Let g e SL(2, Z) and z h- > gz the associated map. It is conformal, one-to-one on EI and maps the real line 
onto itself. The image under such a map of the imaginary axis will therefore be a circle intersecting the real 
axis at right angles. 

Let {a, b} be a pair of coprime integers. Then there are integers p and q such that pa -|- qb = 1 . Therefore 
g = ( b ) € SL(2, Z). The action of g on H maps a point x = iXi, Xt > 0, on the positive imaginary axis 



into the point 



q/1— pq/abx?\ i / 1 
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1/10 



Re[T] 



Figure 4: The function 7tT^({0}) as a function of for — 
subgroup {a, bj c Z x Z whose probability tends to 1 at 



. Each valley of the graph is labeled by the 
= a/b. 



The two parentheses behaves as 1 + 0(t^^) for Ti — > oo. This repeats the statement just made: the image of 
the positive imaginary axis intersects the real line at right angles. The two intersection points are the image 
of T = and oo. Note that, even though the solution p,c|ofpa+c|b = 1is not unique, the form of g was 
chosen so that the image of t = oo and the tangent at this point do not depend of the pair p, q, but only on 
a, b. 

For this particular element g e SL(2, Z), the modular transformation of the probability 7Tt:({1 , 0}) gives 

Because of | [30) the behavior of 7Tt({1,0}) for t — ixt with Tt — > oo fixes the behavior of 7TT^/({a,b}) at 
t' = -g + :^r^. More precisely 



7ra+^J{a,b}) = 1 - ^ Cnq 



with q = e 



(31) 



with the same Cn and as in ijS). Consequently all others 7t({c, d}) with {c, d} ^ ±{a, b) should go to zero 
when T ^ ^. 

Ziff, Lorenz and Kleban fTS\ noticed that the probability n^['L x Z), and therefore 7Tt^{{0}), develop 
oscillations as a function Tr when is close to zero. Their qualitative observation is made quantitative by 
| |3T| . Figure [i] draws the fimction 7tT^({0}) as a function of Tr G [0, j] for ti — (It is sufficient to restrict 
the domain of Tr to [0, as 7tT({0}) — ({0}) and the function f (Tr) — tt^^^ ({0}) is even for a fixed value 
of T|.) The oscillatory behavior is obvious. Each valley of the graph occurs when Tr is a simple fraction a/b, 
for coprime integers a and b and b < 1 and both its width and depth are larger for b smaller, as implied 
by | [3T) . One would see more valleys at a smaller t^. 



4 The limit Q ^ 

The last limit to be taken is not on the geometry, but on the family of models. The partition functions Zq 
on the torus are well-defined for Q e (0,4]. For models with Q in this interval, the Boltzmann weight of 
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any configuration is of the form cQ*/^. The power # is the number I of closed loops for configurations of 
tjrpe {0} and {a, b} and 1 + 2 for those of type Z x Z. As Q goes to zero, the average of the number of loops 
diminishes and configurations with a small number of very long loops are favored. At Q = 0, the set of 
configurations is empty and, consequently, Zq=o and all the partial partition functions Zq=o{{0}), Z.q=o(^x 
Z) and Zq=o({ci, b}) vanish. One may ask what is the homotopy of these very long loops for Q very close 
to zero. Our intuition failed us here. This is why we explored this limit. 

First note that, at Q = 0, the expressions for the partition functions do vanish since, for g — 2 and 



eo = 1, Z({q, b}) — and Z{0} — trivially from ||3j and 1 13 1. This vanishing turns out to be also valid for t 
away from the imaginary axis, but we shall concentrate on the case = for the rest of this section. The 
probabilities 7t(H), H c Z x Z, are therefore the quotient of two quantities that tend to zero when Q — > 0. 
We first expand the partition function Zq (H) around Q = 

Z(H) = fo(H) + efi (H) + yf2(H) + • • • 

where e is a positive number such that 

eo = 1 — e, g=2(l+e) and Q—n^e^. 

As pointed earlier, f o vanishes for every subgroup H. 

The coefficient f i for {a, b} vanishes. Indeed, e appears in ||3| only through g and eo and the first deriva- 
tive with respect to either at e = is easily seen to be zero. The second coefficient fifin, b}) does not vanish. 
The second derivative may be computed by considering the variables g and eo as independent first and 

summing their variations after. Of the three ^— ^r4^^, ^ f ^1'^'^^^ and only the third is not zero. 

Using again Poisson formula, we obtain 

Since Zq (Z x Z) = QZq ({0}) - Tt^e^Zg ({0}) around Q = 0, the probability 7t(Z x Z) can be ignored. The 
computation of Zq ({0}) is shorter as its coefficient f i is non-zero: 

m.m'eZ ^ 

Consequently, 7t({0}) is of order e°, 7t({q, b}) of order e' and 7t(Z x Z] of order e^. At leading order e, they 
are 

7T({0})-1, 7T({Z X Z}) --TT^e^ 

and 

e 2_keZ+l/2l2 aZ+b^T^J^ 



7t({Q,b}) 



Even though loops are very long in typical configurations of models with Q very small, they rarely succeed 
in winding non-trivially around the torus. All sums in 7t({a, b}) are related to elliptic theta functions and a 
compact form is 

7t({Q,b}]~ 



(a2 + b2T2)i 4qe^(q)e3(^) - ^qB^f yq)e3(q) 
where q = qi/(a^+b^T^), 02(q) ^ 2^^o<^^^ q-(-+i> and 03(q) = Lnez q"'- 
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5 Monte Carlo simulations 



Two numerical verifications of the above results were done using Monte Carlo simulations. The first sup- 
ports the claim that Pinson and Arguin's formulae hold for any Q's in the interval (0,4], and not only for 
the integers. The second measures the decay exponent yi /b^ = yi /4 predicted for (1 — 7Tt^({1 , 2})) in the 
limit T= j+iO+. 

Both sets of measurements were done on a family of loop gas models labeled by Q that is known to 
describe the physics of the Q-Potts models when Q is an integer. On figure |5] (a), an Ising configuration 
onaHxV = 4x4 lattice with periodic boundary conditions is drawn obliquely. A (broken) square with 
rounded corners drawn at 45° indicates where the 4x4 lattice is. Because of the tilt, it is useful to label each 
spin by a number 1 to 1 6 to visualize where lie the repeated spins on the boundary. This lattice describes a 
torus with t = i. The basic variables of the loop gas model are the state of the smaller boxes drawn also in 
figure |5] (a). The rectangle with rounded corner that lies horizontally has h x v = 8 x 4 such boxes. 

A Fortuin-Kastele5m (FK) configuration, compatible with the spin configuration, has been chosen in 
figure|5](b). The FK graph is indicated by diagonals in the smaller boxes. The corresponding configuration 
of the loop gas is determined as follows. Note first that two of the vertices of each box is occupied by spins 
of the original lattice. If a FK bond is drawn between them, the state of the box is built out of two quarter- 
circles drawn to avoid the bond. If no FK bond is present, the two quarter-circles are drawn as to prevent a 
bond to appear. Note that the (loop gas) lattice of boxes h x v = (2H) x H has sheared boundary conditions: 
the vertex in the bottom left is repeated in the middle of the top line. This corresponds to t = j + j ■ The 
modulus for the spin lattice (t — i) and that of the loop lattice (t = 2 + 2) ^re distinct, but they lie in the 
same SL(2, Z)-orbit. The Bolztmann distribution on the loop configuration is described in [16J. In [19] a 
simple Metropolis upgrade step is described. The number of steps sufficient for proper thermalization and 
the number of steps between measurements to assure statistical independence are also given there; they 
depend on the model, that is, on Q. 




(a) (b) 



Figure 5: In (a) a 4 x 4 spin lattice with the corresponding loop gas lattice. In (b) an admissible Fortuin- 
Kastele5m configuration with the corresponding loop gas configuration. 



5.1 Models with rational and irrational Q 

We measured the probabilities 7t({0}), n{Z x Z), 7t({1,0}), 7t({0,l}), 7t({l,l}) and 7t({1,-1}), for Q G {1, 

( 2, ir-^^j^)'^]- The four values of Q correspond respectively to percolation, the logarithmic mini- 

mal model CAi(3, 5), the Isrng model and the tricritical Ising model. The two irrational values of Q test our 
claim that l|3]|5| apply to any Q e (0,4]. The cases Q = 1 and 2 were also measured by Argurn 1 1 1 using the 
"spin" models. (The case Q — 1 was first measured in |11|.) We do the measurement here using the corre- 
sponding loop gas models described summarily above. We chose to carry the simulation on a square lattice 
h X h with T = 1. To reduce finite-size effects, the measurements were repeated for h = 4, 8, 1 6, 32, 64, 1 28 
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and the estimates ft(H) were obtained by making a power-law fit of the form 

^(H) -ft^^^(H) = Ci(hx v)^^ 

The results are reported in Table [T] where the 95%-confidence interval is given in the form 0.1 681 14, that is 
0.1681 ± 0.0004. These are statistical errors. The agreement is excellent. Some departure from theoretical 
values is seen for the tricritical Ising model; this was to be expected as this model is closest to the 4-Potts 
model that suffers logarithmic corrections. Results for percolation and Ising agree with 



Model 


VQ 


^({1,0}) 


^({0,1}) 


^({0}) 


7t(ZxZ) 


^({1,1}) 


^({1,-1}) 


CM{3,5) 




0.168114 
0.1680 


0.168214 
0.1680 


0.442715 
0.4429 


0.169114 
0.1692 


0.025812 
0.0258 


0.025712 
0.0258 


percolation 


1 


0.169318 
0.1694 


0.169718 
0.1694 


0.309419 
0.3095 


0.309419 
0.3095 


0.021113 
0.0210 


0.021113 
0.0210 


Ising 


V2 


0.146615 
0.1464 


0.146415 
0.1464 


0.225916 
0.2264 


0.452817 
0.4529 


0.014112 
0.0139 


0.014112 
0.0139 


trie. Ising 


1(V5+1) 


0.130517 
0.1297 


0.130217 
0.1297 


0.196918 
0.1989 


0.5209110 
0.5207 


0.010712 
0.0105 


0.010712 
0.0105 



Table 1: Numerical and theoretical probabilities for six homotopy groups, for the four models correspond- 

ingtoQe{l,(^)^2,(I^)2}. 

5.2 Behavior of 7rT^({l , 2}] close to t = ^ 

We offer only one check of the asymptotic behavior of a 7Tt^(H] on limiting geometries. But it is a non- trivial 
case, 7Tt({1 , 2}), since it probes the exponent yi /b^ obtained in section bI This will be done for the model 
withQ = (l(N/5-1))2. 

The relationship between spin and loop gas lattices described earlier will play here a crucial role. For the 
value of Q under study, there is a loop gas version, but no corresponding spin model. We keep nonetheless 
the name "spin lattice" for the lattice drawn obliquely in figure |5] (a) and use capital letters to give its size. 
(Note that the letter H is also used for the subgroup of the holonomy group. Hopefully this will not cause 
any confusion.) We aim at measuring various probabilities TT-rlH) for spin lattices of size H x V with V = 20 
and H = V/t^ for small ti, that is, for H > V. We choose H = 20, 40, 60, 80, 1 20, 1 60, 320, 640, 960 and 1 280. 
The corresponding value of t = + ix^ are — j and 1 /xt = 1 , 2, 3, 4, 6, 8, 1 6, 32, 48 and 64. To account 
for Ty = J, the bottow row of the spin lattice has to be shifted to the right by C = H/2 sites before being 
glued to the top row. 

For these sizes of spin lattices, the corresponding loop lattices have size h x v (with h and v in small 
letters) given by 

2Vlcm(C + V,H) 2HV 



where 1cm denotes the least common multiple. A shift c, similar to C for the spin lattice, is necessary for 
the loop gas lattice. This shift c is obtained by solving 

cc(-V+C) + |3H = c, a(V+C) + |3H=v 

for a, |3 and c under the constraints a, |3, c G Z and c G [0, h). As an example, the loop lattice h — 1280, v = 
20, c — 620 corresponds to the spin lattice H = 640, V = 20, C — 320. The samples are of 2 x 10^ configura- 
tions for the six smallest lattices and of 1 0^ for the four largest. 
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Figure 6: Numerical values for 7T({i 2}) are plotted, with the theoretical curve and the linear fit obtained 
using the five largest lattices. The size of the dots are larger than the statistical errors. 

The results appear in Table|2] In all cases the agreement for 7t({1 , 2}) is excellent, up to three or four digits. 
This is remarkable considering that one of the lattice linear sizes is very small. Indeed the six largest (loop) 
lattices have v = 20 and finite-size effects should be present. It is also welcome since the measurement of 
Yi requires to take the logarithm of (1 — ft({l,2})]. The slope on figure|6]should be ^^p- at large l/xi. In the 
present case, b — 2 and Yi = ^ = 0.075. Using only the six largest lattices, we extract — 0.0756 ± 0.0005, 
a reasonable agreement, as again the error does not include finite-size effects. 



Appendix 



We prove here the special values | 22p6| of the function C(d, eo). Its definition is 



C(d, eo) = Y_ cos(7Teod2)^L(d/d2) 

dild 

but, if the sum is done over ds = d/d2, it can also be written as 

= ^ cos(7Teod/d3)iJ.(d3). 

dald 

Because the cosine function is even and periodic, the function C(d, eo) satisfies 

C(d, eo) = C(d,-eo) = C(d,eo +2). 
The key relation for proving 1 22p6| is 

C(d, keo) — C(kd, eo) + 6kAd,i C(d, eo), for k prime. 

We first prove it. 



(33) 



15 



1/Ti 


"({0}) 




^({1,0}) 


"({0,1}) 


"({1,1}) 


"({1,-1}) 


"({2,1}) 


"({1,2}) 


1 


0.447717 
0.4480 


0.171016 
0.1711 


0.166316 
0.1661 


0.105115 
0.1047 


0.104515 
0.1047 


0.0024417 
0.00242 


0.0023517 
0.00242 


0.0006514 
0.00059 


2 


0.443417 
0.4429 


0.169015 
0.1692 


0.025913 
0.0258 


0.167916 
0.1680 


0.167616 
0.1680 


0.00010114 
0.00009 


0.000095114 
0.000090 


0.025713 
0.0258 


3 


0.4484110 
0.4486 


0.171018 
0.1714 


0.003912 
0.0039 


0.143017 
0.1424 


0.142717 
0.1424 


0.00000112 
0.000002 


0.000001 12 
0.000002 


0.089216 
0.0894 


4 


0.448517 
0.4480 


0.170916 
0.1711 


0.0006014 
0.00059 


0.104715 
0.1047 


0.104515 
0.1047 


NE 



NE 



0.166016 
0.1661 


6 


0.4107110 
0.4105 


0.156818 
0.1568 


0.00001318 
0.000015 


0.052915 
0.0527 


0.052815 
0.0527 


NE 



NE 




0.3183110 
0.3186 


8 


0.350917 
0.3511 


0.133915 
0.1341 


NE 



0.027013 
0.0268 


0.026713 
0.0268 


NE 



NE 



0.453318 
0.4529 


1 D 


0.151017 
0.1508 


0.057715 
0.0576 


NE 



0.0021919 
0.00220 


0.0021319 
0.00220 


NE 



NE 



0.785018 
0.7850 


32 


0.023113 
0.0232 


0.008912 
0.0088 


NE 



0.00001819 
0.000023 


0.00001718 
0.000023 


NE 



NE 



0.967914 
0.9679 


64 


0.0005015 
0.00053 


0.0002013 
0.00020 


NE 



NE 



NE 



NE 



NE 



0.9993015 
0.99917 



Table 2: Numerical and theoretical probabilities for CM(3, 5). The probabilities '^^[{2-^]) arid 7r(ji _2i) are 
not included because their numerical and theoretical values are < 1 0^. Theoretical values are set to when 
< 1 0^ and NE indicates that no such event was recorded. 

By definition 




Bu summing over all divisors of kd instead of those of d only, the first sum of the second line has added 
terms; the sum Y_ ' is over these spurious terms and restores therefore the equality with the previous line. 
Suppose k is prime. If d has k among its prime factors, all the divisors da of kd that are not divisors of d 
contain k^ as factors. Their Moebius factor \i[d.2) is then and the sum Y.' vanishes. If k is not a prime 
factor of d, then aU di in the sum Y.' are of the form di = kds with da a divisor of d. Then 

LI ( mo kd\ V- /^"^o kd \ 

llT • d. J ^('^^^ = llT • kd3 j ^^^'^^^ 

= -C(d, eo/k). 



We have thus proved 



C(d, eo) = C(kd, eo/k) + 6dAk,i C(d, eo/k). 
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(34) 



Equation | |33) follows if eo is replaced by keo- Both forms are useful. 

The first identity in | 22p6 1 is almost trivial. But it shows how to use | [34) . Suppose d has a repeated 



prime factor, say d = Ic^'. Then k A (kd') = k and the identity l|34) gives C(kd',0) = Clk^d'.O) + 
5kA(kd'),i C(kd',0) — C(d,0). The d's to be studied are therefore those with only distinct prime factors. 
Suppose that d 7^ 1 has I such factors. In the definition of M^fdi), only the number of prime factors is 
important and, if the sum over divisors is replaced by a sum over the number of prime factors in these 
divisors, C(d,0] becomes 

C(d,0) = Y. 1 -Htdi) =(1-1)^=0. 

dild i=0 



Finally C(1 , 0) = m.{1 ) = 1, which proves l|22}. 

The last identity | [26) is the next to be proven. The periodicity of C simplifies its study. If k is an odd 
prime, then C(d, k] = C(d, 1 ). For a given d, choose an odd prime k such that k A d = 1 . Then | [33) gives 

C(d,l) = C(d,k) =C(dk, 1) + C(d,1), 

proving C(kd, 1 ) = 0. This states that 

C(d',l) = (35) 
if d' has a non-repeated odd prime among its prime factors. 



Like above, suppose that d has a repeated odd prime factor, d = k^d'. Then (33 1 and periodicity give 

C(d,l) = C(kd',k) = C(kd',1). 



Removing further factors k if necessary, one can bring these cases back to p5) . The only remaining cases 
are d a power of 2 and d = 1 . If d = 2"^ , n > 1 , then |[33) and periodicity give 



C(2d,1) = C(d,2) = C(d,0)=6d,i = 0. 



A direct calculation gives C(2, 1 ) —1 and C(l , 1 ) — —1, proving p6) . 

Let eo — J now. If d is odd, then all its divisors di will also be and then cos (Tteo dz) — and C ( d, ^ ) =0. 
The periodicity and evenness of C implies also that C(d, |) — C(d, j] for k odd. This allows to use again 
I p3) efficiently. For k odd 

C(d, 1) = C(d, I) = C(kd, i) + 6kAd,i C(d, 1). 

From this point on, the argument is similar to that for C(d, 1 ). The proof of the last two cases {cq — ^ and 
|) uses no new argument and will be omitted. 

These special cases might lead one to think that, for any rational eo, the set {d G N* | C(d, eo) ^ 0} is 
finite. This is false. The cases eo = 0, j , ^ , | , 1 are exceptional in this sense. It is intriguing to note that these 
values of eo are precisely those corresponding to the Potts models with Q =4,3,2,1,0 respectively. (The 
limiting value Q — corresponds to dense polymers.) 
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